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Abstract 

General static solutions for a massless scalar field coupled to a class of effectively 
2-d gravity theories continuously connecting spherically symmetric c?- dimensional 
Einstein gravity {d > 3) and the CGHS model are analytically obtained. They 
include black holes and point scalar charge solutions with naked singularities, and 
are used to give an analytic proof of no-hair theorem. Exact scattering solutions 
in s-wave 4-d Einstein gravity are constructed as a generalization of corresponding 
static solutions. They show the existence of black hole formation threshold for 
square pulse type incoming stress-energy flux, above which trapped surfaces are 
dynamically formed. The relationship between this behavior and the numerically 
studied phase transition in this system |^ is discussed. 



1. Introduction 



Recently, there has been remarkable progress in black hole physics, as we gain valu- 
able insights from model gravity theories such as the CGHS model where the theory 
is exactly integrable [Q]. In these simplified settings, we can address such vexing 
questions as the presumed information loss via Hawking radiation process, black 
hole thermodynamics and related issues in a much more concrete analytical way. 
Even some developments on quantization were manageable 0]. At the same time, 
progress has also been made in 4-d Einstein gravity coupled with various matter 
fields. Especially interesting approach, to name an example, is the study of the 
gravitational collapse in spherically symmetric reduction of 4-d Einstein gravity 
coupled with a single massless scalar field |Q . Several rigorous mathematical results 
and some explicit solutions were obtained so far [0. Moreover, an intriguing 
critical behavior was observed in this system near the onset of the black hole for- 
mation through numerical analysis [|l|. A generic 1-parameter class of solutions S[p] 
decomposes into two phases, depending on the magnitude of p that measures the 
strength of self-gravitational interaction. For p < p*, where p* is the critical value, 
the gravitational collapse of incoming stress-energy fiux is followed by an explosion, 
refiecting back the fiux toward the future infinity. If p > p*, the black holes are dy- 
namically formed and the mass of the asymptotically static black hole Mbh shows 
a universal scaling behavior, Mbh — \p — P*\^ with A ~ 0.37. The analytic expla- 
nation for this phase transition has been only partially successful [Q, although 
in the context of the (1-loop corrected) CGHS model the analytic explanation for 
the similar phase transition with A = 0.5 is available 

As has been well established now, the spherically symmetric 4-d Einstein grav- 
ity and the CGHS model, both of them (effectively) 2- dimensional, can be given 
a unified treatment as a generalized 2-d dilaton gravity theory. To be specific, we 
can construct a class of 2-d gravity theories that continuously connect these two 
interesting theories |TT|. In view of this point, our interest in analyzing these inter- 
connecting theories is at least two- folds; first, we would like to know to exactly what 
extent many interesting results obtained for the CGHS model are the generic lessons 
on gravitational physics and what results are artifacts of the particular choice of pa- 
rameters. Secondly, since the explicit analytic calculations in 4-d Einstein gravity 
are very difficult, we want to find a simpler setting to answer many classical (and 
possibly quantum) questions in general relativity in detail. This can possibly teach 
us more about classical and quantum cosmological questions such as the realistic 
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gravitational collapsing process in general relativity. 

Keeping these considerations in mind, in this paper, we systematically study 
classical general static solutions for a scalar field coupled with various effectively 
2-d gravity theories that continuously connect spherically symmetric 4-d Einstein 
gravity and the CGHS model. The model action we consider here is given as follows 

/ = J d'xy^e''^{R^'^ + ig^''d^<pdf,<p + fie''^ - }-e'^<t>{s~'^ g-^dJOpf) (1) 

and we focus on the case of 5 = 1 for analytic tractability. Here i?*^^-* denotes the 
2-d scalar curvature, (p and / represent a dilaton field and a massless scalar field, 
respectively. 7, A, 5 and are real parameters. The specific choice of these 4 param- 
eters corresponds to a particular gravity theory. Originally, the action introduced 
above was considered as a 2-d target space effective action resulting from string the- 
ory, if we neglect loop corrections. Other interpretations of the action were found 
since then. After {d — 2)-dimensional angular integration, for example, the action 
for spherically symmetric d-dimensional Einstein gravity theories reduces to Eq.(|l|), 
modulo total derivative terms, with 7 = 4((i — 3)/((i — 2), A = 2/{d — 2) and 6=1. 
The dilaton field is related to the geometric radius r of the transversal {d — 2)- 
sphere as r = exp(— 20/(c? — 2)). On the other hand, the CGHS model is recovered if 
we set 7 = 4, A = and 6 = 0. Thus, we see that as long as the pure gravity sector 
is concerned, the CGHS model is exactly the same as c? — 00 limit of d-dimensional 
Einstein gravity. The only difference is the dilaton prefactor in the action for the 
scalar field in the latter theory. 

In section 2, we derive the general static solutions for our model with 6=1 using 
a rather elementary method. Specifically, for general value of 6, we integrate the 
equations of motion, a system of highly non-linear second-order coupled differential 
equations, by finding the relevant number of symmetries of the action, reducing 
the order of differential equations by one. Thus, the result of integration can be 
represented as the conservation of several Noether charges. These equations are 
more analytically tractable and, for 6 = 1 and 2 — A — 7/4>0, i.e. d > 3, they can 
be further integrated exactly to yield closed form expressions. The main novelty here 
is our method of derivation that enables us to get the general static solutions. Unlike 
a well-developed technique of generating solutions for complicated gravity theories 
from vacuum solutions of Einstein equations, our method solves equations of motion 
directly, thereby getting general solutions in a straightforward and simple fashion. 
Our analysis shows that the space of all possible static solutions modulo coordinate 
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transformation is two dimensional (or several sheets of two dimensional space). The 
space is parameterized, roughly, by black hole mass and scalar charge. For non- 
vanishing black hole mass, in 4-d Einstein gravity, we reproduce the exact solutions 
obtained by Janis et. al. |0] long time ago using a solution generating technique from 
the vacuum solutions of Einstein equations. Recalling that solutions we get here are 
general static solutions, our further analysis leads to a new analytic proof of no-hair 
theorem that takes the backreaction on space-time geometry non-perturbatively into 
account. More interesting result is the case of vanishing black hole mass. In this case, 
there is a general relativistic analog of a (non-relativistic) point scalar charge solution 
in the vacuum space-time in each theory. Particularly, these solutions contain naked 
singularities. Although they are solutions for which the total gravity-matter action 
in d-dimensional Einstein gravity 



vanishes, it remains to be seen whether they are stable against classical and quantum 
perturbations. 

In section 3, we find the extension of the static point scalar charge solutions 
into scattering solutions is easy in some cases, namely, the spherically symmetric 
reduction of 4-d Einstein gravity. Utilizing the extension, we construct dynamic 
solutions, which demonstrate that, for the square pulse type incoming stress-energy 
flux, there exists a black hole formation threshold for the incoming stress-energy 
flux Tin- Specifically, in our case, the threshold value of the incoming stress-energy 
flux, T^, that plays the role of p* is found to be 



where G is the gravitational constant and c is the speed of light. If Tj„ > T^, black 
holes are dynamically formed (supercritical case, in the language of phase transition). 
Then, an injection of shock wave type stress-energy flux naturally belongs to the 
supercritical case, for the value of Tj„ is very large during a very short period of time, 
and the corresponding solutions are constructed. Below the threshold, Tj„ < T^, 
the incoming pulse type stress-energy flux gets reflected forward into the future null 
infinity (subcritical case). In supercritical case, we compute the critical exponent for 
our solutions taking the apparent mass of the resulting black hole immediately after 
the formation as an order parameter, getting Ma — \p — p*\^'^- We further indicate 
that our analysis here supports the aforementioned numerical study ]I[ by showing 
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a plausible mechanism that generates a non-trivial scaling relation between Ma and 
Mbh- The difficulties of the similar extension in other generic gravity theories are 
then discussed. Based on the qualitative arguments, we later conjecture that the 
scaling relation for (i-dimensional Einstein gravity for > 4 is Mbh — \p — 
unlike the 4-dimensional (and 3-dimensional 0) results. This dependence of the 
critical exponent on the space-time dimensionality is very similar to the well-known 
cases in condensed matter physics, such as the ferromagnetic transition. 

In section 4, we discuss our results and their possible physical implications in 
various physical contexts. 



2. General Static Solutions 

We start by deriving general static solutions in conformal gauge. The solutions 
we get in this section can also be obtained under a different gauge choice as sketched 
in Appendix. The calculations in Appendix show that, in case of spherically sym- 
metric 4-d Einstein gravity coupled with a massless scalar field, the solutions for 
non-vanishing black hole mass are identical to those found in [0] using a solution 
generating technique. 

2.1. Derivation of General Static Solutions 

The equations of motion we have to solve to get the static solutions are obtained 
from our action by varying it with respect to the metric tensor, dilaton field and the 
massless scalar field; 




l-A 



(2) 




{Dff = 



(3) 



(4) 
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where we define Q = e"^*^ and D denotes the covariant derivative. Since the only 
matter couphng in our case is a massless scalar field, we choose to work in a conformal 
gauge. Thus, we write the metric gais = — exp(2p + ^(j))dx^dx~ . Our convention, 
for the calculational simplicity, is the negative signature for space-like coordinates 
and the positive signature for a time-like coordinate. In spherically symmetric 4-d 
Einstein gravity, for example, this implies /x = —2 using Gauss-Bonnet theorem. 
Under this choice of coordinates, the original action gets simplified to become 

1 = j dx+dx-{md+d-p + ^e^p^l^-^-^'^ + n^d+fd^f), (5) 

modulo total derivative terms. In our choice of the conformal factor, we included a 
contribution from the dilaton field. This contribution, up to a total derivative term 
we threw away, was so chosen to cancel the kinetic energy term for the dilaton field, 
rendering a simplified form of the action. Eq.(^ should be supplemented by gauge 
constraints resulting from the choice of the conformal gauge. They are calculated 
to be 

dlQ - 2d±pd±Q + ^n\d±ff = 0, (6) 
from Eq.(|^) for g±± components of the metric tensor. 

The general static solutions for the equations of motion from the action @ un- 
der the gauge constraints can be found as follows; we can consistently reduce the 
partial differential equations (PDE's) into the coupled second order ordinary differ- 
ential equations (ODE's) by assuming all functions depends on a single space-like 
coordinate x = x^x~ Q Although it is not, in general, possible to have a consistent 
reduction of PDE's to ODE's assuming an arbitrary one-variable dependence, our 
choice turns out to be consistent. This procedure yields the following ODE's 

xh + tl + ^e^P^^-^-^/^p = 0, (7) 
+ P + f (1 - A - 7/4)^ + ^xn'-'r = 0, 



dx 



along with the gauge constraint 



n-2pn + ^n^f'^ 



*This coordinate becomes space- like due the signature choice made here 
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where the dot represents taking a derivative with respect to x. The complete general 
solutions of the above ODE's are the same as the general static Gravity-Scalar 
solutions under a particular choice of the conformal coordinates. The solutions 
under a different choice of conformal coordinates can be obtained from the solutions 
of Eqs-d^) by proper conformal transformations of x^. If we set / = 0, i.e., if 
we consider the pure gravity sector, Birkhoff's theorem ensures that the solutions 
obtained in this fashion are indeed general solutions of the original PDE's. The 
equations of motion other than the gauge constraint can be summarized by an 
action 

1 = j dx{xtlp - t^e^PQ^-^-y/^ - ^xQ^h. (8) 

One can straightforwardly verify that by varying this action with respect to Q, p 
and /, we recover Eqs.(]^). 

We observe that the action has two obvious rigid continuous symmetries. 
First symmetry is clear, for / field appears only through its first derivative. Thus, 
we see that + a is a symmetry for an arbitrary constant a. The second 

symmetry is the transformation x —>■ xe" and p —>■ p — a/2. This symmetry rep- 
resents the scaling invariance of the action (|^). The existence of this symmetry 
is necessary since the original action is invariant under the (local) conformal 
transformations. Since there are three functions p, Q and / that we should solve 
in terms of x, we can integrate Eqs.(|^ once to reduce them to first order ODE's if 
we can find one additional rigid continuous symmetry of the action. The remaining 
symmetry turns out to be x p—>-p — {2 — X — 7/4) ln(l + a)/2 — alnx/2, 

+ a), and / /(I + a)^^~^^^'^. This transformation changes the action 
by a total derivative. We can deduce the form of this symmetry by the following 
physical consideration; the asymptotically fiat spatial coordinate at spatial infinity 
is related to the spatial coordinate that is fiat near the black hole horizon by the 
conformal transformation x^ In familiar from the definition of the tor- 

toise coordinate. The physical system in our consideration is also scale invariant in 
this asymptotically flat region. Our third symmetry is this logarithmic conformal 
coordinate transformation followed by rescaling of fl and / fields, representing this 
asymptotic scale invariance. 

Given these three symmetries, we can construct the corresponding Noether 
charges. 

/o = xQ'f (9) 
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Co = x^ph + \xtl - -yn'p + ^xe^fn^-"--^/^ (10) 
2 4 8 

M + 2co In x = 2xpfi - (2 - A - 7/4)xfi + + ^^xl^^/ (H) 

We can rewrite Eqs. (0) in a form that represents the conservation of Noether charges 
/o, Co and M. By integrating this form of equations of motion, we get the equa- 
tions shown above where /o, cq and M are constants of integration. Additionally, the 
gauge constraint reduces to a condition cq = 0. In the absence of matter fields, mod- 
ulo coordinate transformations, the general solutions of gravity theories considered 
here are parameterized by a single parameter, i.e., the black hole mass. Similarly, 
in our case the gauge constraint kills a redundant degree of freedom, co, for gravity 
sector. As will be clear in the following sections, two remaining Noether charges /q 
and M can be related to a scalar charge and a black hole mass, respectively. 

We can further solve the above equations and this process can be straightfor- 
wardly carried out for 5 = 1. In this case, we can solve p from Eq.(|Tll) and / from 
Eq.(D to find 



P 

and 



/ = /„/f + A. (13) 

where we define y = \q.x and g = l — A — 7/4. Here po and /i are additional con- 
stants of integration. For spherically symmetric reduction of d-dimensional Einstein 
gravity, we have 1 + g = {d — ?>) / {d — 2). Plugging Eqs.(|T2|) and (|TB|) into Eq. (pUj) , 
we get a decoupled equation for fi, an integro-differential equation. 

Our goal is to solve this equation to get as explicitly as possible the closed form 
expression of VL in terms of y. By differentiating it with respect to y we have 

^^a^(M±a^ii^,)(^nr--\fi)l, (15) 

a second order ODE. Integrating Eq.(|T5|) once is immediate; we find 



d^ ^+ d^ fa- 

k9}+'^ = - —\^T^\h. + (16) 
dy dy 

- F (^^) 
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where we define a function h+ and two numbers 



= 20Tqj^y^' + + T M) > 

are introduced. We note tliat A; is a constant of integration introduced when we 
go from Eg. (p!5|) to Eg. (p!6|) . Under the assumption 1 + g > 0, the asymptotic 
analysis of Eg.(|l6D shows that the resulting space-time is asymptotically fiat. In 
this situation, a rather complicated intro-differential eguation ( p!^ can be replaced 
by much simpler Eg. (16). Further asymptotic analysis shows that the eguivalence 
is insured if we identify 

_ / U/2 

4(l + g)^ 

If we properly choose the range of f2 and ^ , we can define the inverse of the function 
Fh_,h+- Then, the integration of Eg.(^) is trivially performed to yield 



where xq is a constant of integration. 

Egs.(|I2]), (|I^) and (|T^ are general static solutions of our problem. Under our 
gauge choice, our solutions are parameterized by 5 parameters, M, /o, po, fi and 
Xq, since the gauge constraint mandates Cq = 0. Among these, /i is just an addi- 
tion of a constant term to /, which is trivial. The constants of integration po and 
Xq represent the degree of freedom in the choice of coordinate system, namely, the 
global scale choice and the reference time choice, respectively. Thus, modulo coor- 
dinate transformations and the trivial /i part, we find the general static solutions 
of Gravity-Scalar action in our consideration are parameterized by two parameters 
/o, the scalar charge, and M, which will be shown to be related to the black hole 
mass. 



2.2. No- hair Theorem and Point Scalar Charges 

In our further consideration of general static solutions, we focus on the cases 

when M > and /o > 0. The restriction on M is motivated by the fact that we 

want the black hole mass to be positive semi-definite. [] The results for negative 

values of /o can be trivially obtained, as is clear from Egs.(0), (0) and ([161) . 

^^In fact, all the solution with M < contains unphysical naked singularities as one can convince 
oneself from the general static solutions. 
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The derivation in the previous section shows that getting the inverse function of 
Fh-,h+ is necessary to determine the relationship between VL and y = \n.x (note that 
we have x = x^x~ in our choice of conformal coordinates). Under our restrictions, 
there are four distinctive cases we have to consider as shown in Figs. la - Id for the 
shape of the function Fh_^h+- In each figure, the physical region to take the inverse 
of Fh^,h+ is shown; since the value of Fh_^h+ relates to some power of geometric 
radius of transversal sphere, via Eg. (p!6D , we require it to vary from zero to infinity. 
Furthermore, we require ^ to be positive for large values of Q, thereby, Fh_^h^. In 
other words, Fh_,h+ = is a natural space-time boundary because it represents a 
space-time point where the transversal sphere collapses to a point. Similarly the 
value Fh_^h+ = corresponds to asymptotic spatial infinities. In every figure, 
Fh-,h+ grows linearly for large, positive showing the spatial infinity regions are 
asymptotically flat. Thus, our choice of physical region is tantamount to considering 
only static solutions with asymptotically flat space-time region at spatial infinities 
and with the conventional choice for the orientation of mappings between Q and 
conformal coordinates there. By making this choice of physical region, we are no 
longer considering parts of general static solutions that have compact range of geo- 
metric radius. The existence of them is clear from Figs.lc and Id. These solutions 
with no asymptotic infinities can be potentially important in some physical settings 
such as cosmological considerations. In this note, however, since we will eventually 
be interested in constructing scattering solutions with asymptotic in/out-regions, 
we restrict our attention to space-time geometries with fiat asymptotic infinities. 



Fig.la. The plot of Fh_^h+ as a function of for h-\- = h- = 0. The region to take the 
inverse of Fh_,h+ is depicted by a bold line. 



10 



lW(l+q) 



-M/(l+q) 

Fig. lb. The plot of F^^^h^ as a function of ^ for /i+ = and /i_ = M/(l + g), a black 
hole geometry with mass M. The region to take the inverse of a bold line, contains 

the point ^ = 0. 

In Fig. la, we have M = and /o = 0, implying h± = 0. Then the scalar field / 
vanishes identically and this is the case that corresponds to the vacuum solution in 
each theory. For non-vanishing value of q, we have 

n={-qk{y-yo))-'^'. (18) 

For d-dimensional Einstein gravity, where we have q = —l/{d — 2), this solution 
represents the expression for the geometric radius in term of conformal coordinates. 
For example, in 4-d case, the above expression becomes = ((Inx"*" + lnx~)/2)^ 
while the examination of Eq. ( [T^ ) shows that In are asymptotically flat conformal 
coordinates. For g = 0, we find 

Inn = k{y-yo), (19) 

a familiar linear dilaton vacuum in pure gravity sector of the CGHS model, as can 
be easily seen after a conformal transform (x^)'^ x^. We can either get this 
expression from Eq.(|l6D or from Eq . ([T8|) taking q ^ limit after replacing ?/o with 
2/0 + !/(#). 

The black hole solutions are recovered if M > and /o = 0, the case in Fig. lb. 
In this case, we have /i+ = and h_ = M/(l + q) and the scalar field / vanishes 
identically. The existence of the apparent horizon, which is the space-time point 
that satisfies d^fl = and, furthermore, is the same as the global event horizon 
in static analysis, is clear from Fig. lb. The indicated physical region to define the 
inverse of Fh_^h+ includes a point ^ = 0. The relation between fl and ?/ = Inx is 
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calculated to be 

ail 



k9}+i - M/(l + q) 
from Eg . (|T7|) , which reduces to 



dy (20) 



2/ 2M, , 

— (r H — — m r — 1 ) = w — wn 

in 4-d Einstein gravity (where g = — 1/2 and = r^) and 

in the CGHS model (where g = 0), reproducing well-known results. (This result 
should be carefully compared with the literature due to our somewhat unconven- 
tional signature choice.) 




Fig. Id. The plot of Fh_^h+ as a function of ^ for a typical /i_ > > case. Dotted 
lines denote a black hole geometry with mass M. The bold line shows the region to take 
the inverse of Fh_^h+- 



In fig.lc, we have M = and /o > 0, resulting h± = fo/ ^2{1 + q). The solutions 
in this case are given by 



/oV(2 + 2q) + 



j dy (21) 



1 ./2(1 + q)k^9?+^<i + fi - fo 

f= , , ) + /i (22) 

/2(l + g) ./2(l + g)W+2, + j2 + j^ 



for n and 



p = i(l + g)lnfi-^y + po (23) 

for / and p, respectively. The integral in Eq.(|21|) can straightforwardly be repre- 
sented in terms of incomplete beta function or hypergeometric function. Examining 
the integral for large f2, we find the asymptotic space-time is the flat vacuum in each 
theory. Apart from the relation between VL and y = \n.x, Eq.(^) is identical to that 
in vacuum solutions in each theory. Thus, the metric gap is fiat in 2-dimensional 
sense and all the information regarding the the curvature of space-time is contained 
in Eg . (pi]) . From Eq.(^^, we see that / diverges logarithmically near f2 = and, 
for large fi, asymptotically falls off like l/Q}^'^. This asymptotic behavior, in d- 
dimensional Einstein gravity, translates to fall-off like 1 / r'^~^ in terms of geometric 
radius r. In between these two limits, / is a monotonically increasing function of 
Vt. The asymptotic solutions for large Vt can also be reproduced directly by taking 
G — s> limit, where G is the gravitational constant. This limit is the same as the 
weak field approximation taking a vacuum metric in each theory as a fixed back- 
ground metric. The wave equation for / in s-wave sector under this fixed background 
geometry has general static solutions of the form / = constant /Vl^^'^ + fi, a solution 
describing a point scalar charge sitting at the origin. Thus, Eqs.(pI|)-(^) are the 
relativistic generalization of the point scalar charge solution. The strong self-energy 
of the scalar field near the origin backreacts to the space-time geometry and, in 
turn, this change softens the power-like divergence of / near the origin into milder 
logarithmic singularity. We also note that the total d-dimensional gravity-matter 
action in d-dimensional Einstein gravity vanishes for these solutions . The physical 
region depicted in Fig.lc does not include the point ^ = 0, a position where the 
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horizon could have been formed. This lack of trapped region implies the singularity 
at Q = is naked. Obviously, black holes are not present in these solutions. 

In Fig. Id, we have generic cases of M > and /o > 0, resulting /i_ > /i+ > 0. 
The field / diverges logarithmically near Q = and the space-time at large Q is 
asymptotically flat. Actually, Fh_^h+ quite rapidly approaches the Fh_^h+ for a black 
hole with mass M as ^ gets large. (See Appendix for the explicit form of / in this 
case.) The space-time geometry is also rather similar to the case of point scalar 
charges, Fig.lc; there are two values of for which Fh_^h^ vanishes. The horizon 
^ = is not included in the physical region to take the inverse of -Fh_,/i+) leading 
inevitably into naked singularities. As a result, black holes are not present in these 
solutions, although M does not vanish. For a better physical understanding of this 
situation, let us consider the case when M » /q > 0, i.e., /i_ ~ + q) and 

a very small positive number proportional to /q. Then, we might try to solve 
the wave equation for / in fixed black hole geometry with mass M, neglecting the 
backreaction of the self-energy of / on space-time geometry. It is well known that 
the static solutions for / in this weak field approximation in curved space-time 
diverge logarithmically at the black hole horizon, ^ = in Fig. Id (for example, 
/ ~ In |r — 2M| in 4-d Einstein gravity near the horizon r = 2M). Since the 
gravitational self-energy of / correspondingly diverges near the black hole horizon, 
the weak field approximation is not valid there. What our calculation shows instead 
is this large self-energy near the horizon backreacts to space-time to cut off the black 
hole from our view and produces a naked singularity in front of the potential black 
hole horizon. Indeed for small but finite /q, the width of the cusp near ^ = h+ 
in Fig. Id gets very small but finite. Only for /o = 0, the cusp disappears and we 
recover the black hole geometry. In other words, when black holes try to carry some 
scalar charge, they end up shutting themselves off from our view, leaving a naked 
singularity. 

To summarize, other than some cosmological solutions with compact range of 
the geometric radius, all solutions with asymptotically flat space-time have naked 
singularities, except for black hole solutions (and vacuum solutions) with identically 
vanishing scalar field /. Thus, we completed the proof of no-hair theorem for all 
model theories in our consideration. Additionally we showed there are relativistic 
analog of point scalar charges in each theory. 
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3. Some Dynamic Solutions 



The static point scalar charge solutions can, in some cases, be simply extended 
to scattering solutions. The s-wave sector of 4-d general relativity is such an ex- 
ample, as shown in this section. The solutions constructed in this fashion describe 
the situations where either black holes are dynamically formed as a result of the 
gravitational collapse or the incoming stress-energy flux gets reflected toward the 
future (null) infinity. 



3.1. The Case of 4-d Einstein Gravity in s-wave Sector 

The spherically symmetric reduction of 4-d Einstein gravity coupled with a mass- 
less scalar field / is the case when 7 = 2,// = —2, A = 1 and 5 = and, as our 
action shows, the gravitational constant G satisfies IQ-kG = 1. Then, the static 
solutions for M = 0, i.e. point scalar charge solutions, are calculated to be 

n = ^{e'^^{Hx/xo)r-4f^), (24) 

1 1 
p= -Infi - -Inx + po, 



„ , , /e^^°^ + /o^-/o , , , , , e^^°ln(x/xo)-2/o , , ^ 
/gSpof] + /2 + e^po ln{x/xo) + 2/o 



where po, fi and xq are arbitrary constants. For simplicity, we take xq = 1 and 
Po = for further discussions. Then, as Q — * oo, the behavior of / asymptot- 
ically approaches to — 2/o/r + fi where the geometric radius r is defined to be 
In this limit, we find r (Ina;"^ + lna;~)/2 and the 4-d metric becomes 
ds"^ — > —dx~^dx~/{x~^x~) — r'^{d9'^ + sin^ 9 dip^) . After a conformal transforma- 
tion Inx^ — > x^, we find that the asymptotic space-time is a fiat Minkowskian. 
We note that under weak field approximation taking a Minkowskian metric as a 
fixed background metric, s-wave sector of static solutions for / is given by / = 
—2/o/r + constant. The physical property of this solution becomes more trans- 
parent if we use r coordinate to describe the geometry of longitudinal space-time 

*We note here that if we couple a massless fermion field and describe the fcrmion s-wave sector 
via a bosonization procedure, the only difference from the above specification of parameters is the 
value of S = 0. 
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instead of conformal coordinates. Then, from Eq.(|2J), we find that the metric is 
given by 

ds^ = Irr^dr^ - r^dO^ + sin" ed(j)^). (25) 

1 + /(f/H 

The above-mentioned asymptotic behavior is clear from this equation. 

The M = case solutions are of particular importance as we can generalize 
the static solutions to dynamic ones. We note the structure of s-wave sector of 
4-d Einstein gravity is very simple as far as the limit G — > is concerned. Then, 
we have a fixed Minkowski space-time and the linear wave equation for / has the 
general dynamic solution of the form / = (/_|_(a;"'") + /_(a;~))/r. Comparing it to the 
static solutions in the same limit, we find that replacing the constant /o in static 
solutions with an arbitrary chiral field generates the general dynamic solutions. This 
consideration suggests that the scalar charge /o, even in general relativistic cases, 
may be a chiral field instead of being a strict constant, similar to what happens in 
G ^ limit. In the framework of general relativity, however, this simple extension is 
not possible in general, since there can be a non-trivial corrections to Eqs.(0) of the 
order of d±fo. In case of asymptotically steady incoming and outgoing stress-energy 
flux, though, the corrections are easily found. Forgetting about global boundary 
conditions, the result of this extension is 

n = ^{Inxf - ^(A;+lna;+ - k-lnx' + go)^ (26) 

111 
p = - InQ - -Inx + - ln(l + k+kJ), 



^ _ j^,lnx - (/i;+lnx+ - A;_ Inx 

In X + (A;+ In x+ — k_ In + go) ' 
where k± and go are constants. We can straightforwardly verify that these so- 
lutions satisfy field equations derived from Eq.(||) and the corresponding gauge 
constraints. Apart from the additional constant term in the expression for p, 
the correction term originating from d±f , Eqs.(p^) are the same as Eqs. (p^) with 
/o = (^o + k-hix~^ — /c+lnx~)/2. The correction turns out to be rather simple 
in this case where the charge /o has terms only up to linear terms in Inx^, which 
are asymptotically flat conformal coordinates near the past or future infinity. The 
asymptotic stress-energy tensor averaged over the transversal sphere in a conformal 
coordinate system f = Inx"'' ( -u = Inx^ ) that becomes asymptotically fiat near 
the past (future) infinity is calculated to be T^^ = k"]^ {Tuu = k'^J)- Thus, k\. are 
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interpreted to be an incoming and an outgoing energy flux, respectively. The con- 
stant go represents the background component of the scalar charge. The qo = 
case of Eqs. (|26|) was reported in mathematics literature ^ (named scale invariant 



solutions) and in physics literature |TB[ where the solutions were used to investigate 
the violation of the cosmic censorship hypothesis. For scattering situations where 
incoming and outgoing flux can coexist, a slightly generalized version proves 
to be useful. Additionally, the presence of go term enables us to consider the time 
evolution of the (multiple) square-type incoming energy pulses by successively glu- 
ing our solutions, unless a black hole is formed in an intermediate stage. To name 
a few other applications possible with our result, we can construct various scatter- 
ing solutions, cosmological solutions and point particle solutions with time-varying 
charge at the origin, depending on the boundary conditions and initial conditions. 

To illustrate a black hole formation in this system, we consider this physical 
situation; in the asymptotic (null) past, we turn on the constant incoming stress- 
energy flux, inject it in a spherically symmetric fashion for a time duration and turn 
it off. After a dynamical evolution, we will be interested in what comes out in the 
asymptotic (null) future. The physical region of space-time in our consideration is 
specified by the requirement i7 > 0, since the angular coordinates should not have 
time-like signature. Thus, the natural boundary is the origin, f2 = 0. In the limit 
G 0, with a fixed Minkowskian background, the situation is represented by the 
solutions depicted in Fig. 2 and given by 

/ = I (27) 

y —kv JJ 

/ = -2k III 

y _ -k{vo+u) JY 

/ = V 
f = ^ VI 

where we use fiat coordinates v = lnx~^ and u = liax~ , thereby Q = = ((u+t>)/2)^. 
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Fig.2. The Penrose diagram shows the space-time region in fixed Minkowskian background. 
The bold fine denotes the space-time boundary = 0. The constant incoming stress- 
energy flux is turned on during the null time v = and v = vq. 

The initial data on the past null infinity were so chosen to describe the turn- 
on of the constant incoming stress-energy fiux r„„ = there at v — and the 
subsequent turn-off at v — vq > 0. We also imposed a boundary condition by 
requiring the field / to be finite along Q = 0. This approximate solution in a 
fixed Minkowskian background has a property that any unbounded amounts of the 
incoming stress-energy fiux are totally refiected off from the origin into the future 
null infinity. However, this picture qualitatively changes as we consider the exact 
solutions focusing on the space-time geometry change due to the stress-energy of 
the scalar field. Under the same choice of the initial data on the past null infinity 
and the boundary condition along Q = 0, we can construct the following relativistic 
solutions. For Q field, we find 

n=\{u + vf I (28) 

n = ^(u + vf - \ev^ II 

O = - k^)v + + k^vof - Ik^Vo + J^f IV 

n^\{{i-k'')v + j^ + k^vof V 



and 



/ = I (29) 



18 



/ 



' v-\-u—kv ' 



v+u+kv > 



/ 



1 / {l-k^)v+u/{l-k'^)+k'^vo-k{vo+u/{l-k'^) ■ 
(i^k'2)v+u/(l~k2)+k'2vo+k{vo+u/{l-k'2) ' 

f = o 



II 
III 

IV 
V 



for / field. The function p can be read off from Eqs.(p6D. We note the limit vq — >■ oo, 
in which we do not turn off the incoming flux, leaves us with the regions I, II and III. 
Our solutions in this limit were previously discovered and discussed by some authors 
[|12| p . The solution in the region VI is rather complicated and, in generic cases, 
can not be obtained using Eqs.(EBD or the conformally transformed version of them 



for the reasons explained later. However, we can obtain the region IV solution 
by matching the III-IV boundary and from a reasonable physical assumption that 
there is no infalling stress-energy flux in the region IV since the asymptotic incoming 



stress-energy flux is turned-off for v > vq. Eqs. (|29|) reduce, after taking the limit 
G ^ (or A; ^ 1) and under the leading order approximation, to Eqs.(0) and the 
geometry is almost Minkowskian, as expected. 




Fig. 3. The Penrose diagram for a typical subcritical case. The boundary 17 = 0, the bold 
line shown above, remains time-like throughout the entire space-time. 

The region I, bounded by the past null infinity, v = and the origin u = —v, 
represents the Minkowski space before the turn-on of the constant incoming flux. 
Regardless of vq, the qualitative properties of our solutions beyond the region I are 
distinctively different for /c < 1 and A; > 1 (notice T„„ = k^). 
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For A; < 1, depicted in Fig.3, all the regions I- VI exist. The region II, bounded 
hy V — Vq, V — and the past null infinity, represents the propagation of the 
incoming particles before any of them hits the boundary $7 = 0. The region III 
represents the region where incoming particles and outgoing particles reflected off 
from f2 = coexist. The net energy flux cancels as a result, and we have a flat 
space-time throughout the region and / fleld, a constant. However, with respect to 
the flat Minkowskian in the region I, the flat space-time in the region 111 is Lorentz 
boosted with relative speed v/c = A;^/(2 — k'^). The region IV contains the outgoing 
particles further propagating toward the future null infinity. In the asymptotic out 
region, which includes the asymptotically flat future null infinity in the region IV 
and the strictly flat region V, there is additional Lorentz boost with respect to the 
flat region III with relative speed v/c — k'^/{2 — k"^) and a translation k'^v^. Due to 
these Lorentz boosts, the original energy pulse looks time-contracted and amplified 
after the refiection off the origin from the point of view of asymptotic in-observers. 
On the other hand, from the point of view of asymptotic out-observers who use 
U {u) = u/{l — k"^) and V{v) = (1 — /c^)v, the outgoing stress-energy flux is Tuu = k"^ 
and the duration of the pulse is fo, exactly the same as what an asymptotic in- 
observer sent in. Thus, the total outgoing energy carried by the outgoing particles 
in the region IV is the same as the total incoming energy thrown in initially. From 
the conservation of energy, we can infer that no further out-going energy flux comes 
out through the future null infinity in the region VI. We also note that the amount 
of the translation of coordinates from the region I to the regions IV and V, k'^v^, is 
proportional to the total injected energy. In summary, for /c < 1, the path of the 
origin remains time-like throughout the whole space-time and incoming particles are 
scattered forward into the future null infinity, just as in the subcritical regime of the 
numerical simulations. 
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Fig. 4. The Penrose diagram for a typical supercritical geometry. The boundary = 0, 
the bold line, becomes space-like for v > and the apparent horizon forms. The details 
shown in the region VI are conjectural. 

If /c > 1, as shown in Fig.4, the path of the origin becomes space-like and form 
a trapped surface in the region II. In this case, the regions III, IV and V disappear 
and our solutions in the region II become exactly the same as the one obtained in 
0. As explained in detail in 0, the resulting space-time is a dynamic black hole 
with increasing mass, for we do not turn off the constant incoming flux up until 
V = vq. Thus, this corresponds to the supercritical phase of this scattering system. 
At /c = 1, which can thus be interpreted as a black hole formation threshold, or 
a phase transition point, the path of the origin becomes light-like in the region II. 
As the numerical studies and the above considerations suggest, the order parameter 
of this system is the black hole mass, which vanishes for A; < 1 and becomes non- 
vanishing for k > 1. Then, the important physical quantity to compute, given our 
exact solutions in supercritical regime, is the critical exponent. The geometric radius 
r = ryi(f ) of the apparent horizon of the dynamic black hole in supercritical case, 
determined by dyV = 0, is calculated to be 



The 1/2 times the value of this corresponds to the apparent mass Ma of the dynamic 
black hole. The linear dependence on v is understandable as it is the time duration 
between the turn-on of the incoming flux and the reference time v. We also note 




(30) 
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that the angular integrated incoming energy flux is Ank^ = jk'^/G. Defining the 
transition point p* = 1 and p = k, we find that the critical exponent in this case is 
A = 1/2 in a scaling relation Ma — {p ~ P*)^- 

In the numerical study |]l[, the scaling relation is given for the asymptotically 
static black hole mass Mbh, measured in the future (null) infinity, and p. In our 
case, on the other hand, since we do not know the solutions in the region VI, the 
scaling relation between and Mbh is not available. However, it is possible to 
show that, in generic cases, we can not choose static solutions in the region VI (in 
the case of G ^ limit, the solutions in the region VI are indeed static as shown 
in Eqs. (pTj) ) , thereby identifying Mbh and M4. To show this, we recall the general 
static solutions in the region VI are given by 

U{u) + V{v) = 2 / (31) 

from section 2 in conformal gauge. Along II- VI boundary, v = vq, we have to satisfy 

^=\{vo + uf-h%l = {r{u)r (32) 

and 

d,n = ^u + ^{l-e)vo. (33) 

By requiring v coordinate to be C°° along the past null infinity, we find V{v) = v. 
Eq.(^) determines the function U{u) via 

/r=r(u) ^fl^p 
F^^) ^^^^ 



where the function r{u) is given in Eq. (^2]) . After some calculations, Eq.(|33|) becomes 

Fh.Mir) = ^(r-h^vor-h^vi (35) 

and this should determine h±. Unless the scalar charge kvQ goes to zero, Eq.(p5|) 
can not be satisfied for any values of h±. Thus, we have shown the above statement 
for generic cases. The only exceptional cases are when the scalar charge vanishes 
after v = vo- A notable example of such exceptions is the shock wave injections [|TU|. 
In this case, we take f — limit for fixed k'^vo, to make the injected energy finite. 
Then, we see that the region II disappears and kvQ vanishes as fo ^ 0. The exact 
solutions, obtained from the above consideration and shown in Fig. 5, are given by 

n = ^{u + vf (36) 
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for the region I, a flat Minkowskian, and 

r 1 

r + 2m In I 1 1 = - (f/ (u) + v) 

2m 2 

along with U{u) = u + Amhi\u/ Am — 1| and m = in the region VI, the 

Schwarzschild geometry. The scalar field vanishes in the whole region except the 
line V = (shock- wave). Recalling Eq.(pOD and the divergence of k, this is an 
example of extremely supercritical cases and we indeed have M4 = Mbh- 



Fig. 5. The Penrose diagram for a shock wave injection along v = 0. The geometry 
for V > vo is described by Schwarzschild metric and the geometry below that is flat 
Minkowskian. 

In supercritical cases, we consequently infer that the space-time geometry in 
the region VI goes through some transient non-static period as long as the scalar 
charge does not vanish for v > vq. The no-hair theorem supports this result. For 
a given finite incoming energy, the theorem insists the only possible candidates for 
the asymptotic final state geometry are black holes with no scalar charge. Thus the 
role of the transient period is to bleach the static component of the residual scalar 
charge. During this process, the apparent horizon that was initially space-like at the 
turn-off time will settle down to a future null direction, the asymptotic final horizon, 
slightly changing its geometric radius. The numerically obtained Mbh — ip ~ V*)^ 
with A ~ 0.37 is very difficult to calculate, as the scaling relation between Mbh and 
Ma gets complicated through this process. In the subcritical cases, we showed that 
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the energy conservation prohibited the emission of out-going flux in the region VI. 
In contrast, the numerically established scaling behavior for the asymptotic black 
hole mass Mbh implies that some energy flux escapes into the future null infinity 
in the region VI for supercritical (especially near critical) cases. 

3.2. The Case of Generic Gravity Theories 

To apply our method of the extension to other gravity theories, we have to get 
a dynamic solution in G — > limit, taking a vacuum metric as a fixed background. 
For example, in d-dimensional Einstein gravity, the linear wave equation in {t, r) 
coordinates in the limit is given by 

- dlr -r\{d-2){d- 4)1/* = (37) 

where / = j* jy'y'^-'^)!'^ and r is the geometric radius. The vacuum solution for the 
space-time satisfies = r*^"^. We immediately find that, unless = 4, the so- 
lutions are quite complicated. Additionally, the static solutions in each case are 

/ = /o/''"'^^^ + /i where /o and /i are constants. Thus, one can imagine setting 
/o = r^'^^^^^'^f* in our static solutions for /* that produces the constant incom- 
ing (outgoing) flux in the past (future) null inflnity, for the relativistic extension. 
Whether this procedure works or not is not yet clear. 

In case of the CGHS model for pure gravity sector, or o? ^ oo limit, the linear 
wave equation is 

- dlr + j/* = (38) 

where / = exp(— A;x/2)/*. Here the geometric background is the hnear dilaton 
vacuum Q = exp(A;x) and x is the asymptotically flat spatial coordinate. The major 
difficulty in this case is the existence of mass term, causing /* to be non-chiral even 
for the asymptotic infinities. 

One interesting point for theories with d > 4 (including d — oo case) is that V^/ 
vanishes as Q — > oo for non-trivial (non-constant) static solutions /. This implies the 
static scalar charge component of sohitions completely decouples from the dynamic 
solutions; the information about the static component can not be contained in the 
initial data specified along the past null infinities, since it vanishes too fast. The 
d = oo case shows this most clearly. The physical frequency spectrum of /* that 
is massive requires the wave length becomes purely imaginary number if we set the 
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frequency to be zero. In fact, the physically allowed frequency is larger than the 
mass in this theory. Thus, the transient behavior of static scalar charge bleaching 
in 4-d Einstein gravity after the turn-off of the incoming energy flux is expected to 
be absent in theories with d > 4. 



4. Discussions 



Our derivation of the general static solutions in section 2 relies heavily on the 
existence of three rigid symmetries of the action. Among these symmetries, there 
is a static remnant from the underlying conformal invariance, which can thus be 
made local. However, the symmetry whose charge relates to a black hole mass does 
not share this property. In the CGHS model with 6 = 0, we have additional rigid 
symmetry Q ^ Q + A where A is a constant. Additionally, this symmetry and 
f ^ f + B can both be made local. These local symmetries resulting from the 
underlying 2-d Poincare current algebraic symmetry are enough to determine the 
general solutions of the CGHS model. One natural question, then, for other gravity 
theories is whether it is possible to find three rigid symmetries that can be gauged 
as what happens in CGHS model. The existence of them will be helpful in finding 
the general solutions in each theory. 

The dynamic solutions in s-wave 4-d Einstein gravity obtained in section 3 are 
interesting from many point of view. First, the dynamics of the point f2 = shows a 
remarkably similar behavior to that of dynamical moving mirror considered in 



for the CGHS model. Trajectories of both points dynamically reacts to the incom- 
ing energy flux. Furthermore, depending on the asymptotic incoming energy flux, 
both show the subcritical behavior where the trajectory is strictly time-like and the 
supercritical behavior where the trajectory becomes space-like. Thus, in some sense, 
the seemingly extraneous introduction of reflecting dynamic boundary (moving mir- 
ror) in the CGHS model makes the model behave more like phenomenologically 
interesting s-wave sector of 4-d Einstein gravity. In the latter theory, the dynamics 
of f2 = is a necessary consequence of the theory, as we showed. Second, it will be 
very interesting to study the quantum field theory on a near critical but subcritical 
geometry. In this case, the incoming and outgoing null coordinates are related by 
very large Lorentz boost, somewhat similar to the infinite red shift near the black 
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hole horizon. Unhke the black hole case, however, no information should be lost 
in the case of subcritical solutions and, additionally, the mode expansion for the 
free scalar field on the asymptotic (null) infinities are well defined. Therefore, the 
near critical and subcritical solutions provide a nice background geometry to study 
gravitational interactions between incoming and outgoing quantum fluctuations. 

Our third interest lies in the phase transition itself. The black hole phase tran- 
sition in the CGHS model is not as difficult as in the spherically symmetric 4-d 
Einstein gravity. The complicated transient behavior is absent in this case, due to 
the simplified dynamics of the model. Therefore, it would be possible to directly glue 
a static dilatonic black hole at v = vq and thereby getting a relation Ma — Mbh if 
we had considered the CGHS model from the the outset. An interesting observation 
in this regard is the critical exponent 0.5 for the scaling relation for the apparent 
mass in s-wave 4-d Einstein gravity is the same as the critical exponent obtained by 
Strominger and Thorlacius in the CGHS model There is a reason for this con- 
nection as suggested in [|lT|. The pure gravity sector of the CGHS model, other than 
being a target space effective action from string theory, can be considered as a lead- 
ing order theory in the l/(i-expansion of the spherically symmetric d-dimensional 
Einstein gravity. We can, therefore, adopt l/(i-expansion and consider the leading 
order behavior in the description of the complex transient process in 4-d Einstein 
gravity. As a zeroth order approximation, we glue a CGHS black hole directly to our 
solutions in the region II to deduce the approximate scaling relation ~ Mbh- 
(See section 3.2.) Thus, the leading order approximation of the exact critical expo- 
nent for square pulse-type incoming energy flux is now calculated to be 0.5. Since 
the next order correction to the critical exponent is expected to be an order of 
1/d = 0.25 and the numerically calculated value is about 0.37, our leading order 
value, 0.5, seems plausible. Furthermore, it is conceivable that the exact critical 
exponent {Mbh — \p ~ P*\^) for pulse-type incoming energy flux in d- dimensional 
spherically symmetric Einstein gravity for d > A would be 0.5, considering the pre- 
sumed lack of transient behavior discussed in section 3.2. In many other cases of 
phase transitions in condensed matter physics, the critical exponent gets the scal- 
ing violation for only lower dimensional cases. It will be an interesting exercise to 
verify this conjecture and, additionally, develop a systematic perturbation theory 
with a dimensionless expansion parameter 1/d to tackle other difficult problems in 
4-dimensional gravity. 
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Appendix. The Alternative Derivation of Static Solutions 

Instead of choosing a conformal gauge, we can choose a gauge used in where the 
general solutions of pure gravity sector in our consideration were obtained to prove 
Birkhoff 's theorem. From now on, we follow the convention of the reference. The 
metric tensor in this gauge is given by 

13^ 

Furthermore, we choose coordinates in such a way that = VL = exp(— 20) and 
require [do, di] = 0. We note that the signature convention, the time component 
being negative, is different from section 2 to closely follow the reference 

The resulting equations of motion can be written as 



dap 



(39) 



— + ^ - ^m^) - ^lidiff + ^(9o/)^] = (40) 



^ + ^ + Y^V{Q) - ^[{djr + g(9o/)^] = (41) 

d,{y'dJ)-do{^n'dof) = (42) 
p a 



along with a gauge constraint 

^ - ln'^of^^f = 0. (43) 

Since we are interested in getting general static solutions, we require dof = and 
dojS = 0. Thus the gauge constraint (|i3| ) is automatically satisfied. Furthermore, 
Eq.(l^) can be directly integrated to yield 

^'^^f = /o^. (44) 
a 



Adding (gg) and (|l]) gives 
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while subtracting them produces 

dilnB = fin^-^^ (46) 

where we define A = a[3 and B = a/ (3. We also set 5 = 1 and ViVt) = ^Q^^. 
Introducing A via A = Q^^^A, we can rewrite above equations as follows. 

4-B = -2ixA (47) 



A = -(2 + 2g + ^) (4^ 



^AA = -(2 + 2g + | 
Ad(j) ^ B 

Here, as in section 2, q is defined as g = 1 — A — 7/4. We can plug Eq.(|47|) into 
Eq.(^) and straightforwardly integrate it once to obtain 

" = ^^^^^ + + ^^^^ " ^^^^ 

where M is the constant of integration. After integrating this once more, we finally 
get 

m^+t = \h+ - B\^^^\h^ + 51^^^+^ = Fh_,h+{B) (50) 
where is a constant of integration and 

1 



h± = ^^rT^(V^^' + 2(1 + q)fS T M). 



We note that Eq. ( |50D is exactly the same as Eq. (0) from section 2 if i? is replaced 
by Using Eq. (^7D , A can be determined as 

jjL ail 

Eqs. (50) and (|5ll) are our main results. They implicitly solve the metric in terms 
of Vt. Then / can be directly integrated using Eq. (^) . 

Choosing i? as a spatial coordinate in favor of VL gives more explicit solutions. 
Then, the metric is computed to yield 
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where the function Q in terms of B is given in Eq. (|50D . Using Eqs.(|^) and (^), 
the scalar field / can be explicitly solved as 



/= , ^° ln||^| + /i. (53) 



Setting B = ^ and M = in the above equation reproduces Eq. (|2^) in section 2. 
The solutions of [0] are recovered if we set 1 + g = 1/2 and A = 1, the case of 4-d 
Einstein gravity, as long as M 7^ 0. We see that it is again possible to derive the 
general static solutions. 
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